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FLUXON INJECTION INTO ANNULAR JOSEPHSON JUNCTIONS 

Alexey V. Ustinov 
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RELATED APPLICATIONS 

[0001] This application derives from Provisional Patent Application Serial Number 
60/283,477, filed April 11, 2001 and incorporated herein by reference, and claims priority 
10 therefrom pursuant to one or more of 35 U.S.C. §119, §120, §365. 

BACKGROUND 

1. Field of the Invention 

[0002] This invention relates generally to Josephson junctions and, more particularly, 
to injection of fluxons into annular Josephson junctions, 

15 2. Discussion of Related Art 

[0003] Long Josephson junctions are interesting systems from the perspective of 
providing a workbench for fundamental investigations of a variety of superconducting 
phenomena, as well as having various applications in cryoelectronics. For an elementary 
introduction see SUPERCONDUCTIVITY by Charles P. Poole, Jr., Horacio A. Farach 

20 and Richard J. Creswick (Academic Press, 1995), pp. 442 - 444 and references cited. 
Long Josephson junctions are also useful for studying basic properties of solitary waves 
(solitons). Solitons of the simplest type are topological kinks and are able to propagate. 
A well-known example of such a soliton is the elementary quantum of magnetic flux Oo 
(also called a fluxon, or Josephson vortex) in a long Josephson junction. See, for 

25 example, A. Barone and G. Paterno, PHYSICS AND APPLICATIONS OF THE JOSEPHSON 
Effect (Wiley, N.Y. 1982); A.V. Ustinov, Physica D 123, 315 (1998). A fluxon in a 
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long Josephson junction can be caused to move along the junction by the application of a 
bias current I B flowing across the junction. The resulting motion of such a fluxon gives 
rise to a dc voltage Vd C across the junction, which is proportional to the fluxon's mean 
velocity u. Thus, a measurement of V dc as a function of Ib provides a useful way to gain 
5 information about properties of the fluxons, including the number of fluxons present. 

[0004] Our primary concern herein is with long Josephson junctions which, for 
economy of language, we refer to simply as "junctions" understanding thereby that long 
Josephson junctions are understood. Explicit descriptions of junctions having other 
shapes will be included when necessary for clarity. 

10 [0005] An important property of an annular long Josephson junction results from the 
quantization of magnetic flux in a superconducting ring. The annular junction is a 
topologically closed system such that the number of initially trapped fluxons is conserved 
and new fluxons can be created only in the form of fluxon-antifluxon pairs. See, e.g., A. 
Davidson, B. Dueholm, B. Kryger, and N. F. Pedersen, Phys. Rev. Lett. 55 2059 (1985). 

1 5 Fluxon motion in annular junctions occurs under periodic boundary conditions and 
without any reflections from boundaries, thereby avoiding many mathematical and 
physical complications that occur for fluxon motion in other junction shapes. One source 
of the interest in investigating annular junctions derives from fundamental aspects of the 
Berry phase effect that arises in annular junctions, (see e.g. F. Gaitan, Phys. Rev. B 63, 

20 10451 1-1 (2001); and V. Plerou and F. Gaitan, Phys. Rev. B 63, 104512-1 (2001)). 
Other sources of interest in annular junctions arise from the phenomena of Cherenkov 
radiation by solitons that can be studied therein, (see, for example, E. Goldobin, A. 
Wallraff, N. Thyssen, and A. V. Ustinov, Phys. Rev. B 57, 130 (1998); and A. Wallraff, 
A. V. Ustinov, V. V. Kuring, J. A. Shereshevsky, and N. K. Vdovicheva, Phys. Rev. Lett. 

25 84, 151 (2000)). Applications with a view towards the development of practical devices 
can also be investigated with annular junctions. 

[0006] Ring-shaped annular junctions have also been proposed as microwave sources 
with high stability and very narrow radiation line width (for example, see U.S. Patent No. 
4,1 81,902 to A. C. Scott). Annular junctions with trapped fluxons have also been 
30 suggested as radiation detectors in which they have an advantage of a stable operation 
point at a finite voltage. (See for example, C. Nappi and R. Christiano, Appl. Phys. Lett. 
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70, 1320 (1997); M. P. Lisitskii et al, NucL Instr. and Methods in Phys. Research A 444, 
476 (2000)). More recently, annular junctions of special shapes have been proposed for 
the creation, storage and manipulation of quantum bits ("qubits") in the form of fluxons 
(see A. Wallraff, Y. Koval, M. Levitchev, M. V. Fistul, and A. V. Ustinov, J. Low Temp. 
5 Phys. 118, 543 (2000)); and fluxon ratchets (E. Goldobin, A. Sterk, and D. Koelle, Phys. 
Rev. E 63, 0311 11 (2001), and Carapella, Phys Rev. B 63, 054515 (2001)). Fluxons in 
Josephson transmission lines, which are discrete analogs of long Josephson junctions, 
have been proposed as on-chip clocks by V. Kaplunenko, V. Borzenets, N. Dubash, and 
T. Van Duzer, Appl. Phys. Lett. 71, pp 128-130 (1997), Y. Zhang and D. Gupta, 
10 Supercond. Sci. Technol., 12, pp 769-772 (1999), D. Gupta and Y. Zhang, App. Phys. 
Let. 76, pp. 3819-3821 (2000), and U.S. Patent Application 6,33 1,805, "On-Chip long 
Josephson Junction (LJJ) Clock Technology", to Gupta et al. 

[00071 A significant problem in utilizing fluxon states in annular junctions is 
preparation of the initial state of the system containing a single or a predetermined 

1 5 number of fluxons. For example, in order to realize a state having a single fluxon, a 
single magnetic flux quantum has to be trapped in the junction, i.e., between its 
superconducting electrodes. The only reliable and reproducible technique for trapping 
fluxons in an annular junction that has been previously known and used requires rather 
exotic and complicated apparatus, namely a low temperature scanning electron (or laser) 

20 microscope. See e.g. A. V. Ustinov, T. Doderer, B. Mayer, R. P. Huebener and V. A. 
Oboznov, Europhys Lett. 19, 63 (1992). Other known methods for trapping magnetic 
flux in an annular junction can be used while cooling the sample below the critical 
temperature of its superconducting electrode(s). These other methods are based on either 
sending a current through an additional specially designed coil placed on top of the 

25 annular junction (see I. V. Vernik, V. A. Oboznov and A. V. Ustinov, Phys. Lett, A 168, 
319 (1992)), or applying a small bias current directly through the junction (see A. V. 
Ustinov, Pis'ma Zh. Eksp. Teor. Fiz. 64, 178 (1996) [Sov. Phys. JETP Lett. 64, 191 
(1996)]; and I. V. Vernik, S. Keil, N. Thyssen, T. Doderer, A. V. Ustinov, H. Kohlstedt, 
and R. P. Huebener, J. Appl. Phys. 81, 1335 (1997)). Unfortunately, the latter techniques 

30 are not sufficiently reproducible and require heating of a junction to high temperature. 
Moreover, fluxons trapped in such ways often suffer from parasitic pinning due to 
Abrikosov vortices which become trapped in superconductive electrodes. Thus, there is a 
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need for a system to inject a single fluxon, or a known number of fluxons, in a controlled 
manner into an annular Josephson junction. The present invention is directed to 
providing such a system. 

5 SUMMARY 

[0008] The present invention relates to a fluxon injection system including injection 
electrodes separated by a distance D in contact with one terminal of an annular Josephson 
junction. Fluxons are trapped on the annular Josephson junction when an injection 
current of sufficient magnitude is injected through the injection electrodes. 

10 [0009] Application of an injection current causes current to flow from one of the 
injection electrodes into a superconducting electrode and across the Josephson barrier. 
The current is collected by another injection electrode on the same superconducting 
electrode so that the total current across the Josephson barrier remains zero. A magnetic 
flux is thus created in the region between the injection electrodes. As the magnitude of 

15 the magnetic flux created by the injection current increases and becomes larger than the 
elementary quantum of magnetic flux, Oo, it may become energetically favorable for a 
compensating negative flux to be created. If the induced flux exceeds 3> 0 , the remaining 
positive flux on the annular Josephson junction can exist in solitary form and become a 
fluxon. The induced flux on the annular Josephson junction is removed when the 

20 injection current is removed. In this case, the compensating negative flux annihilates the 
remaining solitary positive flux and the junction is then free of fluxons. Thus, control of 
the properties of the system, including the current flow between the injection electrodes 
and the electrode spacing, results in controlled insertion of fluxons into the annular 
junction. Such controlled fluxon initial states on annular Josephson junctions can be used 

25 in connection with clock references, radiation detectors, and fluxon oscillators, among 
other applications. Shaped junctions can advantageously be employed along with the 
injection systems of the present invention for initialization of qubits for quantum 
computing, among other applications. 



30 



[0010] These and other embodiments are further described below with respect to the 
following figures. 

-4- 



856727 vl 
M-U584 US 

BRIEF DESCRIPTION OF THE FIGURES 

[0011] Figure 1(a) shows a plan view of an annular Josephson junction with a trapped 
fluxon (dimensions are not to scale). 

[0012] Figure 1(b) illustrates fluxon insertion according to some embodiments of the 
present invention (dimensions are not to scale). 

[0013] Figure 2 shows a schematic view (dimensions are not to scale) of an 
experimentally studied annular Josephson junction having local current injection leads 
according to some embodiments of the present invention. 

[0014] Figures 3(a) and 3(b) show experimentally measured current- voltage 
characteristics of annular junction at injection currents II = 3.37 mA and II = 6.68 raA, 
respectively. 

[0015] Figure 4 shows the measured dependence of critical current of the annular 
junction I c on injection current I L for some embodiments of the present invention. 

[0016] Figure 5 shows the measured dependence of the critical current of the annular 
junction Ic on the applied magnetic field H generated by a current Ih delivered to a coil 
for I H = 0 mA (solid line, a) and for I H = 3.5 1 mA (line with dots, b). 

[001 7] Figures 6a and 6b show numerically calculated dependencies of the 
normalized critical current of the annular junction y c on the injection current amplitude e 
for normalized spacing between injection electrodes, d = 2 and d = 1, respectively, along 
with insets showing the injection current profile. 

[0018] Figures 7a and 7b show numerically calculated current voltage characteristics 
of annular junction normalized injection electrode spacing d = 2 and d = 1, respectively, 
between the injectors and various injection currents. 

[0019] Figure 8 shows the spatial-temporal evolution of the instantaneous normalized 
voltage in the simulated annular system for s=8, d-1, and y=0.4. 

[0020] Figure 9a shows a flux injection system according to the present invention 
utilizing a heart-shaped junction as for qubit initialization, among other purposes. 
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[0021] Figure 9b shows a flux injection system according to the present invention 
utilizing a compact, symmetric and uniformly biased annular junction as for a soliton 
oscillator or clock, among other purposes. 

[0022] Figures 10(a) and 10(b) depict experimental current-voltage plots for two 
5 different values of I L . 10(a): I L =5.27mA. 10(b): I L -10.06 mA. 

[0023] Figure 1 1 depicts the effect of local injection current on the critical current. 

[0024] Figure 12 depicts the dependence of the critical current on the applied 
magnetic field (applied by means of a coil current Ih) for two different values of local 
current I L . 

10 

DETAILED DESCRIPTION OF THE INVENTION 

[0025] The present invention relates to fluxon injection systems. Further, the fluxon 
injection systems according to some embodiments of the present invention are 
demonstrated both experimentally and numerically. Experimental demonstrations of the 

15 fluxon injection systems according to some embodiments of the present invention show 
injection and removal of a desired number of fluxons into or out of an annular Josephson 
junction. Further, a theoretical model of the fluxon injection systems is described. 
Numerical simulations based on the proposed model describing the fluxon injection 
systems show good agreement with the experimental data and provide further incite into 

20 the fluxon injection process and fluxon interaction with small pinning potential remaining 
in the injection region. Several embodiments of the fluxon injection systems according 
to the present invention are described. 

[0026] An annular Josephson junction with a trapped fluxon is shown schematically 
in Figure 1(a) comprising superconducting electrodes la and lb, Josephson tunneling 
25 barrier, 2 and fluxon, 3. The depiction of Figure 1(a) corresponds to the special case in 
which the fluxon, 3, is trapped between the electrodes la and lb when the junction is in 
the process of cooling down through the critical temperature T c of the superconducting 
electrodes. Usually, however, when cooling a Josephson junction through the 
superconducting critical temperature of the electrodes, no such fluxon becomes trapped in 
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the junction. Flux on injection must be accomplished during cooling through T c since it is 
usually not possible to inject a fluxon into the junction when the junction is in its 
superconducting state, that is, below T c . 

[0027] The present invention relates to a fluxon injection system including local 
5 injection of current I L into the superconducting electrodes of the junction by means of 
injection leads attached thereto, schematically shown as 4a and 4b in Figure 1(b). Unlike 
the special case depicted in Figure 1(a), it is assumed for the schematic depiction of 
Figure 1 (b) that there is no magnetic flux trapped in the junction in the absence of 
injection current, that is when Il=0. When the current I L is turned on, current flows from 

10 the injection lead 4a into the superconducting electrode la as 5, and also flows across the 
Josephson barrier, 2. This current II is collected by lead 4b, such that the total current 
across the Josephson barrier 2 remains equal to zero. The current I L generates a local 
magnetic flux O l which, without loss of generality, can be taken to satisfy O l > 0 in the 
region between the injection leads 4a and 4b. As the injection current I L is increased, O l 

15 also increases. As Ol becomes larger than <X>o, it may become energetically favorable to 
have O l compensated by a negative magnetic flux -Oo (with a magnetic field component 
directed outside the ring, see Figure 1(b)). Note that the energy contained in a magnetic 
field increases as the square of magnetic field. Thus, two separated flux quanta O 0 + O 0 
will have lower energy (~ 2O 0 2 ) than a single 2O 0 fluxon (~ 4O 0 2 ). Since the total 

20 magnetic flux in the junction barrier has to remain zero because of the quantization of 
magnetic flux, there also has to appear a positive magnetic fluxon O 0 somewhere else in 
the junction. As soon as the current II gets large enough such that O l > Oo, the induced 
magnetic flux -fOo may exist in the long junction in the solitary form, i.e. as a free fluxon. 
Moreover, if O 0 > 2O 0 , two free fluxons should appear, and so on. 

25 [0028] One anticipates that the existing N free fluxons will interact with the 

remaining magnetic flux (O l - NO 0 ) near the injection leads. An important feature of a 
fluxon injection system according to some embodiments of the present invention is that 
the interaction of free fluxons with the remaining magnetic field (which leads to pinning 
of free fluxons) can be minimized by appropriate choices for the distance D between the 

30 injection leads and for the injection current value II- 



[0029] Experiments have been performed with Nb/Al-A10 x /Nb Josephson annular 
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junction with the mean diameter 2R = 100 \im and the width W = 4 {am. (jxm==10" 6 
meter), both as depicted in Figure 1. The so-called Lyngby geometry furnished with 
additional local current injection leads 4a and 4b as shown in Figure 2 has been used in 
these experiments. The Lyngby geometry is often chosen for experiments as it provides 
5 rather uniform bias current distribution over the junction (Ib in Figure 2) and, at the same 
time, it has proven to be suitable for fluxon trapping at T c . 

[0030] Fluxon insertion using local current injection pursuant to the present invention 
is not limited to the Lyngby geometry of Figure 2, but can be done for essentially any 
annular junction geometry, as discussed further below. 

10 [0031] The junction used for the experimental examples presented herein has a 
critical current density of about 1.1 kA/cm 2 which corresponds to a Josephson length 
Xj » 1 1 jiim and a plasma frequency co p /(27i) « 135 GHz. This implies the ratio 
ItcR/Xj = £^28 of the junction's length to the fluxon's size and the junction width 

W/Xj < 1, i.e., the junction can be regarded as long and quasi-one-dimensional. The local 
15 current I L was applied via 2 (im-wide leads spaced by a distance D = 22 ^m. A magnetic 
field H was applied in the plane of the tunnel barrier (shown as 6 in Figure 2) using a coil 
with a conversion ratio about 0.35 Oe/mA. The measurements were done at the 
temperature of 4.2 K. 

[0032] Measurements were performed in the junction state with no fluxons trapped in 
20 the junction barrier during cooling down. The critical current Ic at I L =0 and H=0 was 
about 7.8 mA. Increasing I L , causes the critical current to decrease and fluxon steps 
appear in the current- voltage (I B -Vdc) characteristics (i.e. bias current vs dc voltage). 
Figures 3a and 3b show two examples of I B -Vd C curves obtained at two different values 
of the injection current II, which was kept on and constant during the course of every 
25 measurement. At Il=3.37 mA (Figure 3(a)), there is a clear single fluxon step with an 
asymptotic voltage of about 64 ^iV. This I B -Vd C curve indicates the fluxon depinning 
current I db to be rather small, about 13% of the junction's critical current Ic measured at 
II = 0. From the top of the step, the junction switches to superconducting gap voltage. 

[0033] Figure 3(b) shows the Ib -Vd c curve obtained at injection current I L = 6.68 mA. 
30 Here we find the double-fluxon step with an asymptotic voltage of about 128 jaV. Again, 
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from the top of this step the junction switches to the gap voltage. We note also the 
remaining depinning current and the single-fluxon step have substantially reduced current 
amplitude. 

[0034] The effect of the local current injection I L on the critical current of the studied 
5 annular junction is presented in Figure 4. Since this dependence was found to be not 
completely symmetric, all four quadrants of the (I c , I L ) plane are presented. The 
similarity of the observed dependence of I c on I L of Figure 4 to the conventional 
Fraunhofer-like pattern of the critical current of a small Josephson junction on magnetic 
field is striking. It can be noted that the single and double fluxon I B -V dc curves observed 
10 in Figure 3 correspond approximately to the first and second minima of the Fraunhofer 
pattern, respectively. 

[0035] Figure 5 presents experimental data on the critical current dependence on the 
applied magnetic field, H generated by a current I H delivered to a coil, both with an 
injection current (I L = 3.5 1 mA) and with no injection current (Il=0). One notes that the 

15 zero-injection pattern shown in Figure 5 (curve (a)) is slightly asymmetric but, in general, 
looks rather as expected for a long annular junction with no trapped fluxons. In contrast, 
the I L = 3.51 mA curve (b) has a pronounced minimum around zero field. Such behavior 
looks very similar to that of an annular junction with one trapped fluxon, as studied in 
detail by I. V. Vernik, S. Keil, N. Thyssen, T. Doderer, A. V. Ustinov, H. Kohlstedt, and 

20 R. P. Huebener, J. Appl. Phys. 81, 1335 (1997), and A. V. Ustinov, B. A. Malomed, and 
N. Thyssen, Phys. Lett. A 233, 239 (1997). Indeed, this is consistent with the single- 
fluxon I B -V dc curve observed in Figure 3(a) which corresponds to a similar range of the 
injection current II- For some ranges at low values of the field H (between coil current 
values between approximately 2mA and approximately 3mA for both polarities), linear 

25 increase of I c with H is observed. 

[0036] A long quasi-one-dimensional Josephson junction is described by the 
perturbed sine-Gordon equation for the superconducting phase difference cp across the 
junction (see A. Barone and G. Paterno, Physics and Applications of the Josephson 
Effect (Wiley, N.Y. 1982)). 

30 <p u ~<p t( =sm<p-a<p t +y +/(*) Eq. (1) 
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[0037] As used in Eq. 1 , x is the spatial coordinate along the junction and t is time, 
measured, respectively, in units of the Josephson length Xj and inverse plasma frequency 
co" 1 , and where subscripts refer to partial derivatives with respect to the indicated 

variable. The coefficient a accounts for the damping due to quasi-particle tunneling 
5 across the junction, while y and f(x) are the spatially-uniform and spatially- varying bias 
current densities, both normalized to the junction's critical current density j c . 

[0038] In the case of an annular junction with no fluxon trapped, solutions of Eq. (1) 
are subject to periodic boundary conditions, § x (x + £) = $ x (x) and <|>(x + i) = §(x) , where 

i - 2nR/X } is the normalized circumference of the junction. 

10 [0039] In order to model the local current injection described above, the spatially- 
varying bias current term can be taken in the form 

f(x) = e{S(x 0 )-b(x Q +d)], Eq.(2) 

where d = D/Xj is the normalized spacing between the injection leads, s is the injected 
current amplitude and 5(x) is the Dirac delta function. The net current J f(x)dx over the 
15 junction is zero. 

[0040] A related model has been studied almost two decades ago by Aslamazov and 
Gurovich (see L. G. Aslamazov and E. V. Gurovich, Pis'ma Zh. Eksp. Teor. Fiz. 40, 22 
(1984) [Sov. Phys. JETP Lett. 40, 746 (1984)). They considered interaction of fluxons 
with an Abrikosov vortex that is trapped in one of the junction electrodes, with its normal 
20 core parallel to the tunnel barrier. The Abrikosov vortex was modeled by 

f(x) = z8 x (x 0 ), Eq. (3) 

where 5 x (x 0 ) is the spatial derivative of the 5 function at x = Xq. Later, the influence of 
the Abrikosov vortex was considered, in that the approach used in the references L. G. 
Aslamazov and E. V. Gurovich, Pis'ma Zh. Eksp. Teor. Fiz. 40, 22 (1984) [Sov. Phys. 
25 JETP Lett. 40, 746 (1984)] and M. V. Fistul and G. F. Giuliani, Phys. Rev. B 58, 9343 

(1998) assumes the condition s«l . In the case described by Eq. (2), the coefficient £ may 
be, in general, arbitrarily large. 
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[0041] Recently, the same problem of an Abrikosov vortex trapped near a long 
Josephson junction was studied in the theoretical paper M. V. Fisrul and G. F. Giuliani, 
Phys. Rev. B 58, 9343 (1998). The Abrikosov vortex configuration considered by these 
authors is substantially the same as that of Aslamazov and Gurovich discussed above (L. 
5 G. Aslamazov and E. V. Gurovich, Pis'ma Zh. Eksp. Teor. Fiz. 40, 22 (1984) [Sov. Phys. 
JETP Lett. 40, 746 (1984)). The main result obtained in ML V. Fistul and G. F. Giuliani, 
Phys. Rev. B 58, 9343 (1998) is qualitatively similar to the findings obtained here in that 
the locally induced magnetic flux generates a stable state of two fluxons with opposite 
polarity, with one of them (antifluxon) pinned by the Abrikosov vortex and another 
10 (fluxon) moving freely in the rest of the junction. For an Abrikosov vortex the 

magnitude of the induced flux cannot exceed d> 0 . Therefore the free energy of the above 
fluxon-antifluxon state is higher than that of the conventional (fluxon-free) state. In the 
examples presented herein, local flux is induced by the external current I L and, thus, may 
exceed (Do. The lowest-energy state then becomes the dissociated fluxon-antifluxon state. 

1 5 [0042] In order to model the process and correctly interpret the obtained experimental 
results on fluxon injection, numerical simulations were performed by solving the partial- 
differential equation of Eq. (1). In the simulations, each 5-function in Eq, (2) was 
approximated by the more smooth function 



20 such that r£ = s. The numerically injected bias is spread over a distance of about 

Through this paper, results obtained with § = 0.5 are presented, which approximate the 
actual experimental case. It has been checked, however, that taking £ « 1 produces rather 
similar results, indicating that the results are not sensitively dependent on the precise 
value of %. In the numerically calculated current-voltage relationships, characteristics are 

25 shown in normalized limits as y(u), were x> is the normalized average fluxon velocity. 
With this normalization u = 1 corresponds to the asymptotic voltage of the single-fluxon 
step. In order to save the computation time, the simulations are performed with the 
reduced normalized junction length t - 10 and the dissipation coefficient a - 0. 1 . 

[0043] Numerical results, in general, turn out to be very similar to the above # 




Eq. (4) 
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presented experimental data. Two different values for the spacing d between the current 
injecting points, d = 2 and d = 1 , have been chosen for simulation. 

[0044] Figure 6 presents the calculated dependence of the critical current y c on the 
injection current amplitude e. As has been already seen in experiment, this dependence 
5 resembles the conventional Fraunhofer pattern of the critical current of a small Josephson 
junction in a magnetic field. The actual length of the corresponding small junction is 
associated with the distance d between the injecting points. As may be expected, the 
overlap between the lobes is larger for larger d, see, e.g. Figure 6(a) and (b). Due to this 
overlap, the minimum value of the critical (fluxon depinning) current between the lobes 
10 is decreasing with d. 

[0045] The calculated curves for normalized bias current y vs normalized voltage v 
for various injection currents (indicated on the plots) are shown in Figure 7. Both single 
fluxon and double fluxon steps can be clearly recognized. Altogether, these curves look 
very similar to those of Figure 3, numerical data show that the steps on curves of Figure 7 
15 account for free moving fluxons under the action of the uniformly distributed bias 
current. 

[0046] It can be seen in Figure 7 that there is a residual pinning of fluxon(s) due to 
the effect of the local bias injectors. According to Figure 6, this pinning is the smallest at 
the local injection current values which lie between the lobes of the y c (e) curve. Thus, for 
20 a given injector spacing d the residual fluxon pinning can be minimized by choosing an 
appropriate value for the injection current s. The smallest pinning can be achieved when 
the lobes of the y c (e) curve join at y c -0, that is for the cases in which d<l . Thus, the 
optimum spacing D between the injectors is about Xj or less (but advantageously larger 
than the width of the injection leads, 1, in order to avoid edge effects). 

25 [0047] This section of numerical results is concluded by Figure 8. It shows a two 
dimensional gray scale plot of the spatial and temporal evolution of the instantaneous 
normalized voltage <p t (x, t) in the annular junction for e = 8, d = 1 and y = 0.4. The 
moving fluxon is recognized as a solitary wave packet moving with a nearly steady 
velocity across the junction. One can see that there is just a very tiny disturbance of the 

30 fluxon motion arising in the region of local current injection. Obviously, this remaining 

pinning may be important only at low fluxon velocities, where the fluxon' s kinetic energy 
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is of comparable scale with the pinning potential. 

[0048] Additional experimental results have been obtained for a reduced spacing 
between injectors of D = 10p.m for comparison with the numerical results presented 
above. Other experimental conditions for the annular Jospehson junction related to the 
5 composition, geometry, temperature, among others are the same as given above in 

connection with the previous experimental results. D=10|i is approximately the same as 
the numerical results presented in Figures 6b, 7b and 8. (The experimental results 
presented herein have d = DA,j -1.1). 

[0049] Measurements were commenced in the state having no trapped fluxons in the 
10 junction barrier during cooling. The critical current, I c , at I L = 0 and H = 0 is about 10.7 
mA. Figure 10(a) and 10(b) give current- voltage plots for two different values of I L - As 
above, the injection current is kept on and at a constant value during the measurements 
presented in Figure 10, Figure 10(a) clearly shows a single fluxon step with an 
asymptotic voltage value about 64 \iV. Figure 10(a) indicates that the fluxon depinning 
15 current is small, about 6% of the junction's critical current measured for I L =0 (I<j p ~ 0.64 
mA). This value is approximately a factor of 2 smaller than that measured above in 
which D was approximately a factor of 2 larger. 

[0050] Figure 10(b) depicts the current-voltage plot for an injection current I L =1 0.06 
mA, or roughly twice the injection current of Figure 10(a). The double-fluxon step is 
20 apparent with an asymptotic voltage of about 128 p.V. 

[0051] We observe in Figure 10 that the current amplitude of the fluxon steps is 
increased by almost 50% in comparison with the previous measurements for D= 22|im. 
At the same time, the remaining fluxon depinning current is reduced. These observations 
taken together indicate higher quality and uniformity of the fluxon states achieved in the 
25 annular Josephson junction. 

[0052] Figure 1 1 depicts the effect of the local injection current on the critical 
current. The similarity to the conventional Fraunhofer patters noted above is not as 
pronounced in Figure 11. 

[0053] Figure 12 depicts the dependence of critical current on the applied magnetic 

30 field (as applied by means of a coil current, I H ). The pattern for II = 0 is more 
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symmetric than that observed in connection with Figure 5 above and is the type of 
functional dependence typical of a long, uniform annular Josephson junction with no 
trapped fluxons. The curve in Figure 12 for Il= 5.75 mA shows a pronounced minimum 
around zero applied magnetic field, a typical behavior for one trapped fluxon. For small 
5 values of H, a nearly linear increase of Ic with H is observed. 

[0054] In summary, local current injected by an injection system according to the 
present invention modulates the critical current of the junction in a manner which is 
similar to the Fraunhofer pattern. Locally produced magnetic fields act on the junction 
region between the injection points, which reacts by creating a magnetic field pattern of a 
1 0 small Josephson junction. By choosing the spacing D between the injecting current leads 
smaller than Xj, one can tune the injection current I L such that the residual fluxon pinning 
is reduced to a very small level. 

[0055] Figures 9a and 9b illustrate additional examples of annular junctions with a 
fluxon injection system according to the present invention. A heart-shaped junction 

1 5 shown in Figure 9(a) can be used as vortex qubit (see A. Wallraff, Y. Koval, M. 

Levitchev, M. V. Fistul and A. V. Ustinov, J. Low Temp. Phys. 118, 543 (2000)) when 
the fluxon quantum state is a superposition of the two lowest energy states in the (upper) 
lobes of the heart. In this case the junction uniformity is required only for the upper 
junction region shown in Figure 9(a), where the fluxon resides. The bottom part of the 

20 heart can be comfortably used for preparing the initial state of the qubit, including 
application of a fluxon injection system. 

[0056] Applications other than qubits often require a compact, highly symmetric and 
uniformly biased annular junction with a fluxon trapped in its barrier. A possible layout 
of such a soliton oscillator or clock is shown in Figure 9(b). The bias current is to be 

25 injected via uniformly distributed thin film resistors (not shown) attached to the 

superconducting film surrounding the junction. The bias current then flows uniformly 
through the barrier and is collected via a contact to the ground plane in the center of the 
ring. It is known that it is nearly impossible to trap a fluxon in such a symmetric 
structure when cooling down through T c . Therefore, local current injection pursuant to 

30 various embodiments of the present invention can be used in this structure as well as in 
any other, see Figure 9(b). The same principle can be also applied to the design of 
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radiation and/or particle detectors based on annular junctions (C. Nappi and R. 
Christiano, Appl. Phys. Lett. 70, 1320 (1997) and M. P. Lisitskii eL al NucL Instr. and 
Methods in Phys. Research A 444, 476 (2000)). Moreover, the working junction area can 
be substantially increased as local flux on injection pursuant to various embodiments of 
5 the present invention are expected also to work in larger two dimensional annular 
junctions, as are typically required for detectors. 

[0057] The fluxon injection procedures described herein offer the advantage of 
reversibility not present in conventional fluxon injection techniques. Reducing the 
injection current to zero causes the fluxon(s) present in the junction to disappear, 
10 returning the junction to the flux on- free state without the necessity of heating the junction 
above its critical temperature. The junction remains in the fluxon-free state until further 
injection current is applied sufficient to re-create fluxons. Thus, practical and convenient 
resetting of the junction to the state of zero fluxons is accomplished. 

[0058] In conclusion, fluxon injection systems have been described herein and 
15 demonstrated both experimentally and numerically. Such fluxon injection systems allow 
for trapping, retaining and removing any desired number of fluxons to, within or from the 
junction. Fluxon injection systems according to the present invention substantially 
simplify the use of annular junctions as oscillators, radiation detectors, vortex qubits 
among other applications. 

20 [0059] Having described the invention in detail, those skilled in the art will 

appreciate that, given the present disclosure, modifications may be made to the invention 
without departing from the spirit of the inventive concept described herein. Therefore, it 
is not intended that the scope of the invention be limited to the specific embodiments 
illustrated and described. 
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